arXiv:1501.03816v2 [astro-ph.EP] 3 Feb 2015 


Draft version February 4, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


THE FEEDING ZONES OF TERRESTRIAL PLANETS AND INSIGHTS INTO MOON FORMATION 

Nathan A. Kaib^'^ & Nicolas B. Cowan^’^ 

Draft version February 4) 2015 

ABSTRACT 

The final stage of terrestrial planet formation consists of several hundred approximately lunar mass 
bodies accreting into a few terrestrial planets. This final stage is stochastic, making it hard to predict 
which parts of the original planetesimal disk contributed to each of our terrestrial planets. Here we 
present an extensive suite of terrestrial planet formation simulations that allows quantitative analysis 
of this process. Although there is a general correlation between a planet’s location and the initial 
semi-major axes of its constituent planetesimals, we concur with previous studies that Venus, Earth, 
and Mars analogs have overlapping, stochastic feeding zones. We quantify the feeding zone width, Aa, 
as the mass-weighted standard deviation of the initial semi-major axes of the planetary embryos and 
planetesimals that make up the final planet. The size of a planet’s feeding zone in our simulations does 
not correlate with its final mass or semi-major axis, suggesting there is no systematic trend between 
a planet’s mass and its volatile inventory. Instead, we find that the feeding zone of any planet more 
massive than O.IMq is roughly proportional to the radial extent of the initial disk from which it 
formed: Aa « 0.25(ai„ax ~ Omin); where a„iin and ttmax are the inner and outer edge of the initial 
planetesimal disk. These wide stochastic feeding zones have significant consequences for the origin of 
the Moon, since the canonical scenario predicts the Moon should be primarily composed of material 
from Earth’s last major impactor (Theia), yet its isotopic composition is indistinguishable from Earth. 

In particular, we find that the feeding zones of Theia analogs are significantly more stochastic than 
the planetary analogs. Depending on our assumed initial distribution of oxygen isotopes within the 
planetesimal disk, we find a ~5% or less probability that the Earth and Theia will form with an 
isotopic difference equal to or smaller than the Earth and Moon’s. In fact we predict that every 
planetary mass body should be expected to have a unique isotopic signature. In addition, we find 
paucities of massive Theia analogs and high velocity moon-forming collisions, two recently proposed 
explanations for the Moon’s isotopic composition. Our work suggests that there is still no scenario 
for the Moon’s origin that explains its isotopic composition with a high probability event. 


1. INTRODUCTION 

An outstanding question in planetary science is the de¬ 
gree to which a planet’s volatile inventory and isotopic 
composition is related to its size and location. In other 
words, how deterministic is terrestrial planet formation? 
In this paper we use ensembles of planet formation simu¬ 
lations to study the statistics of terrestrial planet feeding 
zones. Insofar as the initial planetesimal disks have radi¬ 
ally non-uniform volatile content and isotopic composi¬ 
tion, the width of feeding zones may be used as a proxy 
for the final water inventory and isotopes of planets. 

The water inventory of rocky worlds is i mportant: liq¬ 
uid w ater is the definition of habitability (Kasting et al. 
19931, but planets with too much surface water will have 


no exposed continents, and hence will not benefit from 
the silicate weathering thermostat (Abbot et al. 2012). 
The combined effects of erosion, isostacy, and a substan- 
tial mantle water reservoir mean that planets with water 
mass fractions < 10~^ will not inundate their surface 
(Cowan fc Abbot||2014|. Insofar as a terrestrial planet’s 
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water is delivered in the nor mal course of planet forma¬ 
tion (Raymond et al. 2004), the statistics of planetary 
feeding zones affect the statistics of habitable planets. 

1.1. Theia & the Moon-forming Impact 

Furthermore, understanding the variance in isotopic 
composition predicted from terrestrial planet formation 
is key to deciphering the origin of our Moon. In the 
canonical giant impact hypothesis for the Moon’s ori¬ 
gin, a Mars-mass body (named Theia) strikes the proto- 
Earth in a glancing impact at the tail end of t he gi¬ 


ant impact phase of terrestrial planet formation (Hart- 
mann & Davis||1975l iCameron & Ward| |I976|. 

Llarth, 


This 

impact throws material into orbit around Earth, which 
eventually accretes into the moon we have today. Us¬ 
ing high-resolution smooth-particle hydrodynamic simu¬ 
lations, such a collision has been demonstrated to yield 
a Lunar-mass satellite that is depleted in iron and also 
gives the Earth-Moon system approximately the angula r 


momentum observed today (Canup & Asphaug 2001). 
In such collisions, however, tne moon-forming accretion 
disk around the Earth is largely composed of material 
from the impactor rather than the proto-Earth (Canup] 
2004). This presents a potential problem for the giant 


impact hypothesis because the oxygen isotope compo¬ 
sition of lunar samples has been found to be virtually 


identical to terrestrial rocks (A^’^0< 0.016%o) (|Wiechert 
et al.|20()T). In contrast, meteorites from Mars and Vesta 


are distinctly different from the Earth, with A^'O val- 
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ues of 0.32%o and -0.28%o, re spectively (Franchi et al. 
1999 Clayton & Mayeda||1996 ). More recent work hints 
at a slight ditterence between the Moon and the Earth 
of 0.012%o, but the fact remains that the Moon 

is much more isotopically similar to Earth tha n samples 
from any other large obiect in the solar system (Herwartz 


Because the formation of the Moon requires a relatively 
low-velocity impact between Theia and the proto-Earth, 
it had previously been argued that Theia likely resided 
in an orbit very similar to the proto-Earth and may have 
formed from a simila r region of the planetesimal disk 
( Wiechert et al.||2001 ). If Theia and Earth formed from 
the same pool or plaiietesimals, they may have a similar 
composition, while Mars and Vesta have distinctly dif¬ 
ferent oxygen signatures, owing to their different feeding 
zones in the planetesimal disk. However, the late stages 
of terrestrial planet formation are not a local process. 
Numerical studies of the final giant impact phase of ter¬ 
restrial planet formation have conclusively demonstrated 
that the feeding zones of terrestrial planets are quite large 
and that most bo dies undergo substantial and stochastic 
radial migration (|Chambers| 20011 |Raymond et al.||2004 

Ha 


20051 lO’Brien et al.|| 200 (jj |Haymond et al.|| 200 ti[ 2009_ 
Fischer & Ciesla 2014). As a result, there is no reason 
to expect that Earth and its final impactor formed from 
the same zones of the protoplanetary disk. Hence, the 
identical composition of the Moon seems to rely on a 
potentially improbable coincidence. 

In light of this, there have been recent searches for al¬ 
ternative scenarios to explain the similarities between lu¬ 
nar and terrestrial rocks. Instead of the canonical Mars- 


mass impactor, Canup (2012|) explores the outcomes of 
collisions between a ~0.5 ivl 0 Theia and a ~0.5 Mq 
proto-Earth. For certain impact velocities, the impactor 
makes nearly equal contributions to the Earth and the 
moon-forming disk, resulting in an isotopically similar 
Earth and Moon. This scenario tends to give the Earth- 
Moon too much an gular momentum, how ever. To alle¬ 


Cuk & Stewart 


( 2012 ) suggest that 


viate this problem, 
the system’s angular momentum could have decreased by 
up to a factor of 2 when it passed through a solar evec- 
tion resonance; this allows for a higher impact velocity 
between the proto-Earth and Theia as w ell. If the proto- 
Earth had a rotation period of 2-3 hours, Cuk & Stewart 


( 2012 ) show that high velocity impacts can yield a lunar 


composition domin ated by the proto-Earth. Similarly, 


Reufer et al. (2012) finds that a hit-and-run collision at 
a slightly higher impact velocity can diminish Theia’s 
contribution to the Moon’s composition. 

Another hypothesis for the isotopic similarities of the 
Earth a nd Moon is that they mixed im mediately after the 


impact (Pahlevan & Stevenson 


2007|. The energy from 


the moon-forming impact would have vaporized portions 
of the outer Earth and generated a vapor accretion disk 
initially comprised of the impactor. The disk and outer 
Earth could have mixed and isotopically equilibrated in 
the 10^“^ years that the disk was in vapor form. During 
the time required for mixing, however, the outer portion 
of the disk would have already cooled and begun accret¬ 
ing into the Moon. In addition, in this scenario, one 
would expect refractory elements to exhibit stronger iso¬ 
topic differences between the Earth and Moon since they 


would have condensed into solids faster, yet such a tr end 
is not seen (Zhang et al. | 20 I 2 Armytage et al.|[ 20 I 2 ). 
We may not need to abandon the canonical giant ir 


may not need to abandon tire canonical giant im¬ 
pact hypothesis, however. It is generally presumed that 
isotopically similar Theia and Earth are improbable out¬ 
comes of terrestrial planet formation, but the probabil¬ 
ity has not yet been quantified. While terrestrial planet 
formation simulations have recently been used to con¬ 
strain the original orbit, timing, and collision statistics 
for Theia ( QuarlnsfcLissauer 2014 Jacobson et al. |m4| 


Jacobson &: Morbidelli 2014), few Targe statistical stud¬ 
ies of the potential compositions and feeding zones of 
terrestrial planets (and their impactors) have been un - 


dertaken (Izidoro et al. 2013 Fischer & Ciesla 2014). 


Because the feeding zoiies of terrestriai planets can be 
so large, it is intuitively plausible that Earth acquired 
the di sk-averaged vaMe of due to the planet’s high 


mass (Ozima et al. 
Mars-mass bodies li 


2007|). If this were the case, then 
ce Theia would likely possess similar 


A^^O values, with Mars and Vesta simply being counter¬ 
examples; the unaltered canonical giant impact hypoth¬ 
esis could still be a viable explanation for the Moon’s 
origin. 

With this in mind, we have revisited the problem of late 
stage terrestrial planet formation. We have performed 
150 different simulations modeling the formation of our 
terrestrial planets under different initial conditions. We 
use our simulations to make a statistical study of the re¬ 
lationship between a planet’s final properties—mass and 
semimajor axis—and its accretion history and feeding 
zone. Such relationships can then help us discern how 
often, if ever, the accretion in our simulations results in 
isotopically similar Earth and Theia analogs, but still al¬ 
lows for an isotopically distinct Mars. This will enable us 
to estimate the probability that Theia had an Earth-like 
composition and will help assess the plausibility of the 
canonical giant impact hypothesis. 

The large number of simulations we have performed 
also gives us new insights into the general process of the 
giant impact phase of terrestrial planet formation. In 
particular, we can statistically study the relationship be¬ 
tween the final location of a planet and the location of 
its accretionary feeding zone within the original planetes¬ 
imal disk. With 150 simulations, we will also be able to 
better quantify the level of stochasticity of these feed¬ 
ing zones and how this may affect the delivery of distant 
water-rich planetesimals to the terrestrial planets. 

Our work is organized into the following sections: Sec¬ 
tion [ 2 ] describes the numerical methods and initial con¬ 
ditions employed in our simulations. Following this , we 
have five main sections of results. The first (Section |3.I| ) 
presents a general overview of the relationship between 
the final mass and orbit of a terrestrial planet and its 
feeding zone in the initial planetesimal disk; Section [3.2| 
specifically discusses the fe eding zones for analogs of So- 
lar System planets. Then in |3.3[ we employ several differ¬ 
ent initial hypothetical distributions of A^^O in our sim¬ 
ulations, and we determine how often Earth-like Thieas 
are formed alongsi de isotopically distin ct Mars analogs. 
With the work of Canup (2012) and ICuk & Stewart! 


( 2012 ) in mind, we next use our simulations to quan¬ 
tity the mass and imp act velocity dist ributions of Theia 
analogs in SectionIn Section [T5] we then make pre- 
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dictions for the isotopic composition of Venus relative to 
Earth based on the distributions employed in the 

study of Theia’s composition. We summarize the con¬ 
clusions of our work in Section ID 

2. NUMERICAL METHODS 

We have performed three ensembles of 50 simulations 
modeling the giant impact phase of terrestrial planet for¬ 
mation in the solar s ystem with the M ERCURY hybrid 
integration package (Chambers 1999). Each ensemble 
makes different assumptions about the distribution of 
planetesimals as well as the orbits of Jupiter and Saturn. 
The first ensemble, “circular Jupiter and Saturn” (CJS), 
assumes that Jupiter and Saturn begin with nearly cir¬ 
cular orbits (e < 0.01) on their present-day semimajor 
axes. We begin with 100 self-interacting embryos orbit¬ 
ing between 0.5 and 4.0 AU on nearly circular (e < 0.01), 
nearly coplanar (i < 1°) orbits. These embryos’ semima¬ 
jor axes are spaced evenly, and their masses decrease with 
semimajor axis to yield a surface density profile. 

In addition, each simulation includes 1000 smaller plan¬ 
etesimals that do not gravitationally interact with each 
other but do interact with all of the other simulation bod¬ 
ies. These planetesimals are given a fixed mass, and their 
initial semimajor axis distribution is adjusted to yield a 
surface density profile as well. Like the embryos, 
the planetesimals are initially placed on nearly circular, 
coplanar orbits. Arguments of perihelion, longitudes of 
ascending node, and mean anomalies are randomly gen¬ 
erated from uniform distributions for all bodies. Our 
disk mass of ~5 Mq is split evenly between embryos and 
planetesimals. Finally, we set the radii of embryos and 
p lanetesimals by a s sumin g a bulk density of 3 g/cm^. 


Raymond et al. (2009) demonstrated that terrestrial 
planet formation simulations with an initially circular 
Jupiter and Saturn consistently produce planets near 1.5 
AU that are much more massive than the real Mars (Fig¬ 
ure Q. Initializing Jupiter and Saturn on eccentric orbits 
quickly truncates the planetesimal disk via secular res¬ 
onances, resultingin Martian analogs more akin to the 
real Mars (Figure]^. We therefore include a second set of 
simulations, “eccentric Jupiter and Saturn” (EJS). These 
simulations start Jupiter and Saturn with moderately ec¬ 
centric orbits (e = 0.1) at their present-day semimajor 
axes. The distributions of embryos and planetesimals in 
EJS is identical to that employed in CJS. 

Another failure of many terrestrial planet formation 
simulations is that they fail to form the two most mas¬ 
sive planets between 0.7 and 1.0 AU as our own solar 
system has. However, the radial concentration of plane¬ 
tary mass seen in our own solar system can be replicated 
if terrestrial plan et formation begins from a narrow an - 
nulus (Figure [si (jMorishima et al. 2008 Hansen 2009). 
We mimic this m our last ensemble, “annulus" (AiNiN). 
As in CJS, we assume nearly circular coplanar orbits for 
Jupiter and Saturn (e < 0.01, i < 1°). However, in this 
ensemble all of our embryos are confined to initial semi¬ 
major axes in an annulus between 0.7 and 1.0 AU. This 
also crudely replicates the plane t formation scenar io put 
forth in the Grand Tack Model (|Walsh et al.|2011|) . The 
ANN simulations contain 400 ernbryos, each with a mass 
of 0.005 Mq, on nearly circular, coplanar orbits with ran¬ 
dom mean anomalies, longitudes of ascending node, and 
arguments of perihelion; no planetesimals are present. 


This fin al ensemble ha s identical starting conditions to 
those of Hansen (2009). 

We choose these three sets of initial conditions because 
they are similar to those employed in numerous popular 
past studies of terrestrial planet formation, and we wish 
to see how objects’ accretion histories vary under differ¬ 
ent assumptions. Each set has shortcomings and there 
are many other pot ential initial conditions that we do 


not e xplore [e.g. ][](Raymond et al. 2009 Izidoro et al. 
2014). It is beyond the scope of the present work to ad¬ 


vocate one set of initial conditions over another. Rather, 
we employ these different initial conditions to explore the 
possible parameter space of the solar system’s terrestrial 
planet formation. 

For all simulations, we integrate orbits with a time 
step of 6 days for 200 Myrs. Particles are removed from 
simulations if they pass within 0.1 AU of the Sun or if 
they go beyond 100 AU. The computing time to complete 
this large suite of simulations is quite significant, and 
we will make our simulations publicly available at http: 
//nathcinkaib. com, 

3. RESULTS 

3.1. Feeding Zone Trends 

We first consider whether there are general trends 
between the location and width of a planet’s feeding 
zone and its final properties, namely semi-major axis 
and mass. While previous works have studied feeding 
zone stochasticity and correlation wit h final orbit (e.g., 
Chambers||20011 |Raymond et al.|j2004[ |Fischer & Ciesla| 
2014D , our large number of simulations and large parti¬ 
cle number per simulation allow us to search for general 
statistical trends that may have been less obvious in past 
work. 

We pay special attention to planets with attributes 
similar to Solar System terrestrial planets and have de¬ 
fined categories for Venus, Earth and Mars analogs in our 
simulations. Venus analogs have final semimajor axes be¬ 
tween 0.6 and 0.75 AU and masses greater than 0.5 M^. 
Similarly, planets are designated as Earth analogs if 0.8 
< a < 1.2 AU and m > 0.5 Mq, while Mars analogs 
are defined as 1.3 < a < 1.7 AU and m > 0.05 Mq. 
The ranges of a and m employed to define our analogs 
are quite liberal, but this must be done to generate a 
statistically useful sample of analogs of each class (Fig¬ 
ures mg]). 

In uigures and we plot the planetary feed¬ 

ing zones vs. inrm semi-major axes for our three ensem¬ 
bles of simulations. The location of a planet’s feeding 
zone is the mass-weighted mean semi-major axis of its 
constituent planetesimals; the feeding zone width is the 
mass-weighted standard deviation of its constituent plan¬ 
etesimals’ semi-major axes. 

Planets that end up close to their star tend to form 
from material close to the star, and vice versa, but the 
relation is not 1:1 because the initial disk has a finite ex¬ 
tent and the material spreads out. Moreover, individual 
planets have feeding zones that span a large fraction of 
the initial disk, and the width of the feeding zone does 
not exhibit a strong trend with final-semi-major axis. 

Nor does the width of a planet’s feeding zone correlate 
with its mass (Figure]^. In fact, while low-mass planets 
have a variety of feeding zone widths, planets more mas- 
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Fig. 1. — The final mass and semi-major axis for the CJS ensem¬ 
ble (initially circular orbits for Jupiter and Saturn). The colored 
boxes indicate analogs of Venus (green), Earth (blue), and Mars 
(red). 



Fig. 2. — The final mass and semi-major axis for the EJS ensem¬ 
ble (initially eccentric orbits for Jupiter and Saturn). The colored 
boxes indicate analogs of Venus (green). Earth (blue), and Mars 
(red). 

sive than O.IM 0 have approximately the same size feed¬ 
ing zone for each ensemble. This behavior is even more 
striking if one considers the unweighted standard devia¬ 
tion of initial semi-major axes, as shown in Figure!^ In 
other words, planets accrete planetesimals from a fairly 
consistent feeding zone, but the few accreted planetary 
embryos tend to be more of a crapshoot; this has im¬ 
portant consequences when we consider Theia analogs, 
below. 

The constant feeding zone size of terrestrial planets dif¬ 
fers between our three ensembles: Aacjs = 0-7(1) AU, 
AttEjs = 0.36(7) AU, and Aoann = 0.086(5) AU 
(the values for the un-weighted standard deviation are 
Aacjs = 0.68(8) AU, Aoejs = 0.36(5) AU, and 
Aoann = 0.086(5) AU). The radial mix ing statistic 
(|Chambers|2001||Quintana & Lissauer|2014|), Sr, is 0.57, 
0.39, and 0.27 tor the CJS, FJb, and AJNJN ensembles, re¬ 



Fig. 3.— The final mass and semi-major axis for the ANN en¬ 
semble (0.7-1.0 AU annulus of planetary embryos). The colored 
boxes indicate analogs of Venus (green), Earth (blue), and Mars 
(red). 



Fig. 4.— The location (mass-weighted mean) and width (mass- 
weighted standard deviation) of a planet’s feeding zone is plotted 
vs. its final semi-major axis. The diagonal gray line denotes a 1:1 
correspondence between planetesimals and the planets final loca¬ 
tion. 

spectively. 

The feeding zone widths are proportional to the initial 
disk size from which the planets form. The CJS ensem¬ 
ble starts with planetesimals and embryos at 0.5-4.0 AU 
for a width of 3.5 AU. The EJS ensemble effectively has 
an initial annulus of 0.5-2.5 AU (width of 2.0 AU) since 
planetesimals beyond 2.5 AU are rapidly ejected by res¬ 
onances with the eccentric gas giants. Finally, the 0.7- 
1.0 AU annulus configuration has a width of 0.3 AU. 
Given these three ensembles, it therefore appears that 
the width of a planet’s feeding zone is 20-30% the width 
of the initial planetesimal disk. 

3.2. Feeding Zones of Solar System Analogs 

For each analog planet we have gone back and looked 
at the relative contribution that planetesimals and em- 
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Fig. 5.— The location (mass-weighted mean) and width (mass- 
weighted standard deviation) of a planet’s feeding zone is plotted 
vs. its final semi-major axis. The diagonal gray line denotes a 1:1 
correspondence between planetesimals and the planets final loca¬ 
tion. 



'O.O 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

planetary mass [Earth masses] 

Fig. 7.— The mass-weighted standard deviation of initial semi¬ 
major axes is plotted vs. its mass. The symbol shape denotes the 
ensemble of simulations from which the planets were drawn. For 
each ensemble, the median feeding zone width for planets larger 
than 0.1 M 0 is shown as a black line. The colored symbols are 
analogs of Venus (green), Earth (blue), and Mars (red). 



Fig. 6 .— The location (mass-weighted mean) and width (mass- 
weighted standard deviation) of a planet’s feeding zone is plotted 
vs. its final semi-major axis. The diagonal gray line denotes a 1:1 
correspondence between planetesimals and the planets final loca¬ 
tion. 

bryos in our initial disk make toward its construction. 
Weighing all bodies by their mass, we are then able to 
see which regions of the protoplanetary disk are sampled 
most heavily during the construction of Venus, Earth, 
and Mars analogs. The results of this analysis are shown 
in Figure panel A shows the CJS simulations. We see 
that Mars analogs in these simulations have a distinct 
composition relative to Venus and Earth analogs. In the 
previous section we saw that planets forming in a disk 
tend to all have the same size feeding zone, so it stands 
to reason that Earth and Venus (0.3 AU apart) will have 
more overlap than Earth and Mars (0.5 AU apart). Al¬ 
though there is a stochastic nature to terrestrial planet 
formation, this plot shows that Mars analogs consistently 
form from material orbiting further from the Sun. This 
is less true for Venus and Earth analogs: although on 



Fig. 8 .— The (unweighted) standard deviation of initial semi¬ 
major axes is plotted vs. its mass. The symbol shape denotes the 
ensemble of simulations from which the planets were drawn. For 
each ensemble, the median feeding zone width for planets larger 
than O.IM 0 is shown as a black line. The colored symbols are 
analogs of Venus (green), Earth (blue), and Mars (red). 


average Earth analogs form from more distant material 
than Venus analogs, the two planets’ cumulative distri¬ 
bution functions overlap significantly. In addition, we 
find that analogs of all three terrestrial planets typically 
receive comparable contributions of material beyond 2.5 
AU, which is presumably water-rich (|HayashilI981l|Ray-| 

ImoAd et al.||2004|. -*-"-' 

We also searcn our simulations for Theia analogs, 
which we define to be the last body with m > 0.1 Mq 
that struck an Earth analog. (Previous studies have 
found that our moon could be formed with impactors 
less massive than 0.1 M 0 , but it is clear from Figures!^ 
andj^that our simulation’s mass resolution limit strongly 
affects these smaller bodies’ feeding zones.) We see in 
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Fig. 9. — The relative contribution of each region of the plan- 
etesimal disk to the final compositions of Venus, Earth, Mars, and 
Theia analogs are show with green, blue, red, and hashed CDFs, 
respectively. The median CDF of each class of analogs is shown 
with a solid line, and the I-ct range of CDFs seen amongst our 
simulations for each analog class is marked with the shaded region. 
Panels A, B, and C show the results from the CJS, EJS, and ANN 
simulations, respectively. 


Figure]^ that the composition of Theia analogs varies 
much more than that of planetary analogs. In some 
cases Theia analogs are built from material more dis¬ 
tant than Mars analogs, while in other cases they have 
compositions dominated by material nearer to the Sun 
than Venus analogs’ compositions. Part of the reason for 
this variation is that Theia analogs are simply smaller 
bodies on average. Thus, their accretion history sam¬ 
ples fewer regions of the disk, and an ensemble of Theia 
CDF curves displays more variance than planetary CDF 
curves. However, by definition Theia is a planetary mass 
body that is quasi-stable at best. Eventually it becomes 
unstable and strikes the Earth. Theia’s very wide CDF 
swath in Figure suggests that before impact, Theia 
can be quasi-stable interior to the Earth (implying a com¬ 
position dominated by inner disk material) as often as it 
is quasi-stable exterior to the Earth (which would imply 
a greater fraction of outer disk material in its composi¬ 
tion) . This fact also helps explain the large compositional 
variance of Theia analogs. 

The EJS simulations are shown in Figure [^. In these 


simulations, material from outside 2.5 AU makes nearly 
no contribution to the compositions of the planets. This 
is because the eccentricities of these distant planetesi- 
mals are rapidly excited via strong secular and resonant 
perturbations from Jupiter, which quickly leads to their 
ejection before they can be accreted by the growing ter¬ 
restrial planets. In addition, we see that Venus, Earth, 
and Mars analogs all have narrower CDE swaths with 
less overlap, indicating a more deterministic outcome of 
planet formation. However, Theia analogs still display a 
very large compositional variance compared to the plan¬ 
ets, as some Theia analogs are largely formed from mate¬ 
rial beyond Mars whereas others are mostly formed from 
material between Venus and Earth. 

Finally, we investigate the ANN simulations in Figure 
HP- Here we find that Venus and Earth analogs typi¬ 
cally have very similar compositions with material at all 
distances making nearly identical contributions to both 
planets. Eurthermore, there are less dramatic compo¬ 
sitional differences between Earth and Mars analogs in 
this ensemble compared to CJS and EJS. Mars has only 
a slight bias toward more distant material, and there is 
overlap with the CDF swath for Earth analogs. In fact, in 
12.5% of the runs that contain Earth and Mars analogs, 
the mean location of the Martian feeding zone is actually 
closer to the Sun than the Earth analog’s feeding zone’s 
mean location (in stark contrast to our CJS or EJS sim¬ 
ulations). Once again, we see a substantially higher level 
of variance for the composition of Theia analogs com¬ 
pared to the planets. 


3.3. Oxygen Composition of Earth, Theia, and Mars 

Because of the chaotic nature of the final stage of 
planet formation, an Earth-like oxygen isotope composi¬ 
tion for Theia is considered improbable. While two past 
studies of the values within a single planet forma¬ 

tion simulation suggest this is the case, the probability 
that Theia possessed an Earth-like oxygen isotope com¬ 
position has never been estimated, even though terres¬ 
trial p lanet formation simulations have the potential t o 
do so (Ziegler et al. 2006 Pahlevan fc Stevenson||2007 ). 
Given an initial distribution of among embryos 

and planetesimals, it is easy to predict the A^’^0 values 
of the final planets by weighting each accreted body’s 
A^^O contribution by its mass. This approach is chal¬ 
lenging because we have virtually no constraints on the 
functional form of the initial A^^O distribution in the so¬ 
lar nebula. Nevertheless, we can test different functional 
forms and use the difference between Earth’s and Mars’ 
fractionation as constraints on our chosen hypothetical 
distributions. (This approach of employin g hypotheti¬ 
cal distribution s is similar to that taken by Quintana fc:| 
Lissauer (2014|) when studying water delivery during ter¬ 
restrial planet formation.) With the constraints provided 
by Earth’s and Mars’ isotopes, we can then attempt to 
estimate the probability that the final major impactor 
with the Earth (Theia) acquires an Earth-like distribu¬ 
tion of oxygen isotopes. 

Below we investigate several different simple families 
of initial distributions of A^^O on our simulations that 
contain Earth, Mars, and Theia analogs. To increase 
the fraction of analyzed simulations, we further relax our 
criteria for Mars analogs to simply be the first planet 
orbiting beyond an Earth analog (as defined above). This 
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yields 68 different simulations for analysis (37 from CJS, 
16 from EJS and 15 from ANN). 


3.3.1. Linear Distribution 

The simplest non-uniform initial distribution we can 
assign for is a linear function in a. In each simula¬ 

tion that forms Mars, Earth, and Theia analogs, we im¬ 
pose on our initial planetesimals and embryos a distribu¬ 
tion in that varies linearly with a. In each simula¬ 

tion, the slope and intercept of the distribution is chosen 
so that the final A^^O differs by 0.32%o between Mars 
and Earth. Thus, each simulation has a unique initial 
A^^O distribution tuned to replicate the observed differ¬ 
ence between Earth and Mars. The slopes of these dis¬ 
tributions are shown in Figure [Toj^. For our simulations 
with a disk of planetesimals between 0.5 and 4 AU (CJS 
and EJS), the A^'^0 gradient is of order l%o per AU. It is 
much steeper for ANN, since the annulus has a smaller 
spread in initial semimajor axes and forces Earth and 
Mars analogs to have feeding zones that overlap more. 

With these A^^O gradients chosen, we can now look 
at the distribution of A^^O values for the Theia analogs 
in each simulation. These are shown in Figure [To)3. (We 
plot the absolute magnitude of A^^O, since we are more 
concerned with the difference in A^^O than its sign.) As 
can be seen, the difference in |A^^O| for Earth and Theia 
is almost always greater than that observed in the real 
solar system (a ssuming the Moon is composed of mate¬ 
rial from Theia ( CaiiuppOOJ I). The median |A^^O| value 
of Theia analogs m the CJb, EJS, and ANN simulations 
is 0.17, 0.13, and 0.12%o, respectively. Out of the 37 
systems we study in the CJS set, none of them manage 
to generate a Theia analog with |A^’^0| below 0.016%o. 
In the EJS simulations, only 1 of 16 systems manage to 
produce an Earth-like Theia. Similarly, in the ANN sim¬ 
ulations only 1 of 15 systems generate a Theia analog 
with |A^^O| < 0.016%o. The failure to generate Earth¬ 
like Theia analogs in these simulations is surprising. In 
all but two simulations we fail at generating a Theia ana¬ 
log with an Earth-like oxygen composition in spite of a 
smooth linear A^^O designed to produce the observed 
difference between Earth and Mars. This suggests that 
either the canonical moon formation scenario is unlikely, 
or that the initial A^^O distribution in the solar system 
was very different from a linear distribution. 

The canonical collision scenario for the Moon’s origin 
predicts that at least 70% of the M oon’s material will be 
derived from Theia (|Canup||2004|). More recent alterna¬ 
tive collision scenario s successtully decrease that fraction 


to values of ^5-40% (Cuk & Stewart|2012 


irt||201 


Reufer et al. 


20121. From our results in I'igure [TOp, for each 'I'heia 
analog we can easily calculate its maximum contribu¬ 
tion to the Moon’s composition allowed by the measured 
lunar isotopic composition. The distribution of these al¬ 
lowable contributions is shown in Figure |10p . Here we 
see that the large majority of our Theia aiialogs are pro¬ 
hibited from contributing more than ~40% of the Moon’s 
mass. The median value of the allowable Theia contri¬ 
bution is 9.7%, 13%, and 13% for CJS, EJS, and ANN 
simulations, respectively. 

We can also search for a possible trend between |A^^O| 
and the time of the moon-forming impact. For instance, 
in the case of the CJS simulations, we may expect late- 



Fig. 10.— A: Assuming an initial distribution of |A^^O| among 
planetesimals and embryos that varies linearly with heliocentric 
distance, a histogram of the |A^^O| distribution slopes is shown 
for all of our simulations that form Earth, Mars, and Theia 
analogs. Each distribution is tuned so each simulation yields 
|A^^O| = 0.32%o when comparing the final oxygen isotope com¬ 
positions of Earth and Mars analogs. B: Using the distributions 
specified in panel A, the | A^^O| values of Theia analogs in each sim¬ 
ulation are calculated and shown here. The shaded regions marks 
values of |A^^O| below 0.016%o. C: Given the Theia lA^^O] val¬ 
ues in panel B, the maximum allowable contribution to our Moon’s 
mass is calculated for each Theia analog. The CDFs of these allow¬ 
able contributions are shown here. CJS, EJS, and ANN simulation 
results are displayed with the solid, dashed, and dotted CDFs, re¬ 
spectively, in all three panels. 


impacting Theia analogs to come exclusively from the 
distant outer part of the planet-forming disk, since ac¬ 
cretion proceeds much more slowly in this region. If this 
were the case, then late-colliding Theia analogs may show 
especially large differences between themselves and the 
Earth analogs. To search for such correlations, in Figure 
ddwe plot the |A^^O| value of each Theia analog against 
the time it collides with the Earth analog. As can be seen 
in the plot, there is no obvious correlation between im¬ 
pact time and how isotopically Earth-like a Theia analog 
is. 

The results shown in Figure [T0| do not necessarily in¬ 
dicate that Earth and Theia anSlugs have fundamentally 
different feeding zones. Even if the two bodies had the 
same feeding zone we would expect a large fraction of 
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Fig. 11. — The | A^^0| values of Theia analogs are plotted against 
the simulation time at which they collide with the Earth analog. 
Circles, squares, and triangles mark the CJS, EJS, and ANN sim¬ 
ulation sets, respectively. 




Fig. 12. — A: The | A^^O| values of Theia analogs in each simula¬ 
tion are shown. The dashed line shows the CDF for a set of artificial 
Theia analogs built by simply resampling the feeding zone of each 
Earth analog, while the solid line is the CDF for the actual Theia 
analogs constructed via accretion in our simulations. The shaded 
regions marks values of |A^^O| below 0.016%o. B: CDFs of the 
maximum contribution to the moon’s mass that each Theia analog 
can make. Our artificial Theia analogs are shown with the dashed 
line and our actual Theia analogs constructed via accretion in our 
simulations are shown with the solid line. 

our Theia analogs to have different values simply 

because they are smaller and sample the same pool of 
planetesimals fewer times. Thus, the population of Theia 
analogs should have a higher variance in than the 

population of Earth analogs as long as the Theia analogs 
have lower masses. 

We can easily test whether Theia’s lower mass is the 
dominant source of the large A^^O values among Theia 
analogs. To do this, in each simulation we remove the 
Theia particles from our Earth analog, so we have an 
Earth feeding zone that includes only non-Theia parti¬ 
cles. Next we redraw particles from our modihed Earth 


feeding zone until we have drawn the equivalent mass of 
the original Theia analog that we just discarded. With 
this done, we have constructed an artificial Theia ana¬ 
log whose feeding zone is forced to be the same as the 
Earth analog. Using our modified Earth analog and our 
artificially constructed Theia analog, we can then pre¬ 
dict a A^^O for the artificially built Theia analog just as 
we did in Figure [TO] Moreover since we are constructing 
our new Theia analogs without dynamical accretion, we 
repeat this resampling process 100 times for each simu¬ 
lation to build up a smoother A^^O distribution. The 
results of this process are shown in Figure [T2]4. Here we 
compare our artihcially constructed Theia analogs with 
the real Theia analogs that were generated via accre¬ 
tion in our simulations. We see that the Theia analogs 
built from resampling Earth’s feeding zone have signif¬ 
icantly smaller |A^^O| values (median of 0.054%o) than 
our original Theia analogs (median of 0.147%o). In addi¬ 
tion, only 2.9% (2/68) of our original Theia analogs had 
|A^'^0| < 0.016%o, whereas 20% of our artificially built 
Theia analogs have |A^'^0| < 0.016%o. We can also again 
look at the maximum contribution that our set of Theia 
analogs are allowed to make to the moon’s composition 
in Figure [T^. We see that our artificially constructed 
Theia analogs can make up a significantly larger fraction 
of the moon’s mass (median value of 29.4%) compared to 
the original set of Theia analogs generated by our simu¬ 
lations (median value of 10.9%). 

Although Figure [T^ shows that some of the variance in 
the A^’^0 values ofTheia analogs is due to the bodies’ 
smaller masses, it also demonstrates that Theia analogs 
tend to have fundamentally different feeding zones than 
Earth analogs, and the differences in feeding zone loca¬ 
tion between Earth and Theia is what drives most of the 
A^^O difference between the Earth and Theia analogs 
in our simulations. This can be more clearly seen in 
Figures miA -C. In these panels we calculate the differ¬ 
ence between the mean formation semimajor axis of each 
Earth and Theia analog pair. A CDF of these differences 
is plotted for each simulation set. This is done for the ac¬ 
tual Theia analogs built in the simulations as well as the 
artificially generated Theia analogs used to make Fig¬ 
ure [12] If Theia and the Earth were derived from the 
same material the two CDFs would overlay each other 
in each panel. Instead, the real Theia analogs show a 
significantly greater spread in formation semimajor axis. 
The plots indicate that although Theia is roughly equally 
likely to form from material interior to the Earth as it is 
to form from exterior material, its feeding zone is usu¬ 
ally distinct from the rest of Earth. Using a K-S test, 
we compare the artificial Theia CDF with the real Theia 
CDF for each simulation set. This comparison yields p- 
values of 0.13%, 0.13%, and 7.8% for the CJS, EJS, and 
ANN simulations, respectively. 

3.3.2. Step Function 

In contrast to smooth linear A^^O distributions, we 
can also explore initial A^’^0 distributions designed to 
give sharp contrasts between outer and inner terrestrial 
planets. To do this, we model A^^O distributions with a 
step function. To build our step function, we first choose 
a contrast for the inner and outer values of A^^O. Next, 
we set the heliocentric distance at which the step occurs 
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Fig. 13. — The difference between the mean formation semima¬ 
jor axis of each Earth and Theia analog is shown. (Negative Aa 
implies that Theia formed from material that is interior to Earth’s 
material.) The solid line is a CDF of the semimajor axis difference 
between Earth analogs and the Theia analogs formed in our simula¬ 
tions via accretion. The dashed line is a CDF for the Theia analogs 
that are artificially constructed by simply resampling the Earth’s 
feeding zone. The results of the CJS, EJS and ANN simulations 
are shown in panels A, B, and C, respectively. 


in order to produce a difference in |A^^O| for Mars and 
Earth analogs that is as close as possible to 0.32%o. (Be¬ 
cause the smallest masses in our simulations are of order 
10 “^ M 0 , it is not always possible to exactly match a 
0.32%o difference with a step function.) 

The results of this approach will clearly depend on the 
contrast we choose between the inner and outer values 
of A^^O. For this reason, in each simulation we test a 
range of contrast values. Guided by the spread of mean 
A^^O values seen among the chondritic mete orite groups 
( Burbine fc O’Brien|2004 Ozima et al.|2007 ), we explore 
contrasts between U and l 0 %c. in some cases, it is impos¬ 
sible to attain a difference near 0.32%o for Earth and Mars 
analogs’ values no matter what heliocentric distance the 
step occurs. We therefore only consider simulations that 
have attained a Martian |A“0| of 0.32 ± 0.03 relative 
to Earth. In Figure [T4|4, we show the fraction of our 
simulations that have met this criterion as a function of 
the chosen contrast. One can see that a |A^^O| a dif¬ 
ference of 0.32%o between Earth and Mars cannot be 


attained for small step sizes. This is not surprising, as 
no difference between the two planets’ compositions will 
occur if the step contrast is zero. For CJS and EJS sim¬ 
ulations, the optimal contrast is between ~1 and ^4%o. 
Beyond 4%o the number of Earth and Mars analogs with 
|A^^O| near 0.32%o falls off again. The reason for this 
is that while there is often a large amount of overlap be¬ 
tween the feeding zones of Earth and Mars analogs in 
these simulations, some simulations yield fairly distinct 
feeding zones. As a result, it sometimes becomes impos¬ 
sible for the two analog planets to have a small |A^^O| 
value if a large contrast is chosen (no matter what the 
heliocentric distance of the step is). Alternatively, the 
ANN simulations require a larger contrast in the initial 
|A^^O| distribution to attain Earth and Mars analogs 
with |A^^O| near 0.32%o, because planets forming from a 
narrow annulus have nearly entirely overlapping feeding 
zones. 

For step function contrasts in which at least half of our 
simulations match the Martian |A^^O|, we use this same 
step function to calculate the |A^^O| value for Theia 
analogs. Of the simulations that yield a Martian |A^^O| 
of 0.32 ± 0.03, we plot the fraction of Theia analogs that 
have |A^^O| < 0.016%o in Figure [l4|3. From this plot we 
see that the | A^^O| values for Theia analogs are typically 
greater than 0.016%o, regardless of the position and con¬ 
trast of the step. Examining all the versions of the initial 
|A^^O| we employ, the mean fraction of Theia analogs 
with |A^'^0| < 0.016%o is 5.5%, 7.2%, and 4.4% for our 
CJS, EJS, and ANN simulations respectively. These fre¬ 
quencies of Theia analogs with Earth-like oxygen signa¬ 
tures are comparable to those found when we assumed a 
linear distribution of A^^O. Thus, an Earth-like composi¬ 
tion for Theia is an unlikely outcome of terrestrial planet 
formation for either family of initial A^^O distributions. 

Once again, we also determine how much material each 
Theia analog can contribute to the Moon’s composition. 
We calculate the median value of this allowable contri¬ 
bution for each step function contrast we employ for the 
initial A^^O distribution. The results of this are shown in 
Figure 14P. Similar to our linear A^^O distributions, we 


find that typically Theia material must not account for 
more than 5-15% of the Moon’s composition. The excep¬ 
tion to this occurs at small step function contrasts, where 
Theia is sometimes allowed to contribute up to 40% of 
the Moon’s mass. However, as can be seen in Figure T4]4, 
these small step contrasts have difficulty consistentlyre- 
producing Mars’ isotopic difference. 


3.3.3. Random Distribution 

The fact that the A^^O of Vesta is lower than Earth, 
while that of Mars is higher is suggestive that the initial 
distribution of A^^O did not monotonically increase or 
decrease with heliocentric distance. This would invali¬ 
date both linear trends and step functions for the initial 
radial A^^O distribution. The opposite extreme is to as¬ 
sume that there is no correlation between initial A^^O 
values and initial semi-major axis. To perform this ex¬ 
ercise, we generate a A^^O value for each initial particle 
in our simulations by randomly drawing from a Gaus¬ 
sian distribution centered on A^^O = 0. By combining 
the mass-weighted A^^O values of our original particles 
we then calculate the final |A^^O| of Earth, Mars, and 
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Fig. 14. — A: The fraction of Mars analogs that are able to satisfy 
I A^^O| = 0.32±0.03%o (compared to the system’s Earth analog) as 
a function of the chosen jump in | A^^O| in the initial | A^^O| step- 
function distribution (for all simulations containing Earth, Mars, 
and Theia analogs). B: The fraction of systems that are found 
to have Earth and Theia analogs with |A^^O| within 0.016%o of 
each other as a function of the chosen jump in | A^^O| in the initial 
|A^^O| step-function distribution. C: Given the predicted |A^^O|, 
the maximum allowable contribution of Theia to the Moon’s mass 
is calculated. The median value of the maximum allowable Theia 
contribution is plotted as a function of the step size in in the as¬ 
sumed initial |A^^O| distribution. CJS, EJS, and ANN simulation 
results are shown with black circles, blue diamonds, and red stars, 
respectively, in all three panels. 


Fig. 15. — A: The distribution of |A^^O| for Mars analogs if 
we assume that the initial A^^O values of planetesimals are ran¬ 
domly drawn from a Gaussian distribution (for simulations that 
contain Earth, Mars, and Theia analogs). The shaded region 
marks 0.29%o< |A^^O| < 0.32%o. Initial A^^O values of em¬ 
bryos and planetesimals were redrawn 10000 times to build our 
statistics. B: The distribution of |A^^O| values for Theia analogs 
for the same simulations as in panel A. The shaded region marks 
|A^^O| < 0.016%o. C:Given the Theia |A^’^0| values in panel B, 
the maximum allowable contribution to our Moon’s mass is calcu¬ 
lated for each Theia analog. The GDFs of these allowable contri¬ 
butions are shown here. CJS, EJS, and ANN simulation results are 
displayed with the solid, dashed, and dotted GDFs, respectively, 
in all three panels. 


Theia analogs. We rerun this random generation of 
10,000 times on each ensemble (CJS, EJS, and ANN) to 
build up our statistics. For each ensemble, the stan¬ 
dard deviation of our initial A^’^0 values is adjusted so 
that the median absolute difference between the A^^O 
of Earth and Mars analogs is 0.32%o. To meet this re¬ 
quirement, the standard deviations of our initial |A^^O| 
distributions are set to 5.80, 4.86, and 2.82%o for the 
CJS, EJS, and ANN simulation sets, respectively. 

In Figure 15 \, we display the CDF of |A^^O| values 
for our Mars analogs for the 10,000 random oxygen iso¬ 
tope generations in CJS, EJS, and ANN. As can be seen, 
the |A^^O| of Mars analogs varies substantially between 
^0.01 and ^l%o with the median near 0.32%o. In Fig- 
ure|l^, we plot the CDFs of the | A^^O| values for Theia 
analogs in these same simulations. In all three sets of our 


simulations, there is only a tiny probability (~2% in each 
simulation set) that a Theia analog will have a |A^^O| 
below 0.016%o. The level of variance in A^^O required to 
generate the observed oxygen difference between Earth 
and Mars analogs almost always also yields Theia analogs 
that differ substantially from Earth. Thus, once again we 
find it quite unlikely that an Earth-like composition for 
Theia is produced via terrestrial planet accretion. 

Figure [T5| J shows the maximum contribution that each 
Theia analog can make to the Moon’s mass. Again, only 
small mass fractions are typically allowed. The median 
values of the allowed Theia contribution to the Moon’s 
mass are 3.4%, 4.1%, and 4.5% for the CJS, EJS, and 
ANN simulations, respectively. 

In our previous hypothetical A^^O distributions with 
an explicit dependence on heliocentric distance, the fi- 
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nal values of the planets and impactors are greatly 

affected by the level of radial mixing within the proto¬ 
planetary disk during terrestrial planet formation. How¬ 
ever, when one assumes no dependence on heliocentric 
distance (as is the case here) the final compositions of 
planets and their impactors are largely independent of 
the simulation’s dynamics. Instead, the expected devia¬ 
tion from the disk-averaged value for each planet 

and impactor is a function of the number of planetesimals 
and embryos that the body accreted during the simula¬ 
tion; i.e., its mass. This is just a manifestation of the 
Central Limit Theorem, and since we tune our initial 
A^^O to typically yield substantial differences between 
Earth and Mars, it is no surprise that Theia analogs 
(which tend to be less massive bodies) usually have A^^O 
values that are at least as large as the simulations’ Mars 
analogs. 

3.3.4. Bimodal Values 


One last scenario that we explore for the initial A^^O 
values in the protoplanetary disk is one where all em¬ 
bryos have one fixed A^^O value and all planetesimals 
have another. One could imagine such a situation aris¬ 
ing if bodies’ A^^O values are affected by their mean for¬ 
mation time within the solar nebula or different parent 
body processing. In this case, each protoplanet’s A^^O 
value would be set by the fraction of its mass contributed 
by planetesimals vs embryos. (Note that the ANN simu¬ 
lations are excluded from this analysis since they began 
with only embryos.) 

To set the A^^O difference between planetesimals and 
embryos in each simulation, we simply calculate the plan- 
etesimal fractions of the Earth and Mars analogs and 
then fix the A^^O of planetesimals to yield a Martian 
|A^^O| of 0.32%o. The distribution of A^^O differences 
between planetesimals and embryos in our simulations 
is shown in Figure [T6|4. We see that planetesimals typ¬ 
ically must have |ZV^O| values that are ~2%o different 
from the embryos in order to match the observed isotopic 
difference between Earth and Mars. 

With the A^^O of our initial embryos and planetesi¬ 
mals tuned, we can now determine the |A^^O| value pre¬ 
dicted for Theia in each simulation. These are shown in 
Figure pT)3. We see that there are nearly no examples 
of Theia analogs with Earth-like oxygen isotopes, and 
the typical isotopic difference between Earth and Theia 
is 0.1-0.2%c. Thus, if the Moon’s composition is dom¬ 
inated by material from Theia, our simulations predict 
that its observed isotopic composition is a very improb¬ 
able event. 

This fact is emphasized in Figure [T^ where we calcu¬ 
late the maximum mass fraction that each Theia analog 
is allowed to contribute to the Moon’s mass. We find that 
^90% of our Theia analogs are prohibited from consti¬ 
tuting more than ~30% of the Moon’s mass, and the 
median allowable mass fractions are 7.2% and 9.5% for 
the CJS and EJS simulations sets, respectively. 

However, all of our simulations fix th e initial mass ratio 
of planetesim als to embryos to be 1:1. [Jacobson fc Mor-| 


bidelli (2014) investigated planet formation with various 
embryoAo-planetesimal ratios and concluded that the to¬ 
tal mass was likely dominated by embryos at the start 
of the giant impact phase of planet formation. Further¬ 
more, given the shorter collision timescales in the inner 





Fig. 16.— A: Assuming two different initial for planetesi¬ 

mals and embryos, a CDF of the |A^^O| differences between plan¬ 
etesimals and embryos required to reproduce the observed Martian 
isotopic value is shown for all of our simulations that form Earth, 
Mars, and Theia analogs. B: Using the distributions specified in 
panel A, the | A^^O| values of Theia analogs in each simulation are 
calculated and shown here. The shaded regions marks values of 
|A^^O| below 0.016%o. C: Given the Theia |A^^O| values in panel 
B, the maximum allowable contribution to our Moon’s mass is cal¬ 
culated for each Theia analog. The CDFs of these allowable contri¬ 
butions are shown here. CJS, EJS, and ANN simulation results are 
displayed with the solid, dashed, and dotted CDFs, respectively, 
in all three panels. 


protoplanetary disk, it’s likely that the true ratio of plan¬ 
etesimals to embryos varied with distance, and presum¬ 
ably the number of planetesimals would be greatest in the 
outer part of the disk. If this is the case, one could imag¬ 
ine Mars accreting many more planetesimals than the 
Earth and Theia. If so, an isotopically identical Earth 
and Theia may not conflict with an isotopically distinct 
Mars in this case. This could also potenti ally allow for 


an is o topically uniform inn er solar nebula (Javoy et al. 


2010 Dauphas et al. 2014), with Mars being polluted 


by iriore distant planetesimals. However, this scenario is 
untested and should be the subject of future studies. 


3.4. Potential Solutions: Massive Theias and High 
Velocity Impacts 

In the previous section, we found that regardless of the 
initial A^^O we impose on embryos and planetesimals, it 









































12 



7 


Fig. 17. — The distribution of 7 values seen among Earth analogs 
in our simulations (7 being the ratio of Theia’s mass to the com¬ 
bined mass of Theia and the proto-Earth at the time of impact). 
Results from the CJS, EJS, and ANN simulations are shown with 
the black solid, dashed, and dotted CDFs, respectively. The com¬ 
bined CDF for all Earth analogs is marked with the thicker gray 
line. 



Fig. 18. — The cumulative distribution of impact velocities be¬ 
tween Earth and Theia analogs in the ANN simulations. Theia 
analogs are split into three different mass bins: m = 0.025-0.05 
Mq (solid line), m = 0.05—0.1 M 0 (dashed line), and m > 0.1 M 0 
(dotted line). Impact velocity is calculated in terms of the mutual 
escape velocity of the Earth and Theia analogs. 


is very unlikely for Theia analogs to be isotopically simi¬ 
lar to Earth analogs. This presents a significant issue for 
the canonical moon formation scenario involving a Mars- 
mass impactor, since it generates a moon primarily com¬ 
posed of the impactor rather than the Earth. However, 
Canup (2012) suggests that this issue can be resolved 
it 'i'heia had a mass comparable to that of the proto- 
Earth. In this case, both the Earth and moon-forming 
disk are a roughly even mixture of the proto-Earth and 
Theia. (This scenario relies on the angular momentum of 
the Earth-Moon system later decreasing via a n evection 
resonance with the Sun ( Cuk fc Stewart[2012 1 .) 

With our large number of terrestrial planet formation 
simulations, we can estimate the statistical likelihood 
that Theia’s mass was comparable to the proto-Earth. 
To do this, we simply look at the distribution of mass 
ratios for Earth analogs struck by Theia analogs in our 
simulations. This distribution is shown in Figure 17 


where we plot the parameter 7 , which is the ratio o: 
Theia’s mass to the combined mass of Theia and the 
proto-Earth at the time of impact. T o mix the Earth 
and Moon evenly enough, Canup (2012|) finds that Theia 
must have had 7 > 0.4. in EigureTTTT we see that such 
collisions are not found in any of our simulations. Out 
of the 104 Earth analogs generated in our collisions, the 
largest recorded 7 is 0.325, and only 8.7% of our Earth 
analogs experienced impacts with 7 > 0.3. Late im¬ 
pacts with 7 > 0.4 must be exceedingly rare, implying 
that a comparably massed Theia and prot o-Earth is a 
very unlike l y even t. This result agrees with Jacobson fc| 
Morbidelli (2014), who also found that major mergers 
between protoplanets with similar masses are rare. 

Instead of comparable masses fo r The ia and the proto- 
Earth, both 


Reufer et al. 


( 2012 ) and 


Cuk & Stewart 


(2012) invoke a higher impact velocity between 'I'heia 


and the proto-Earth to produce an Earth-like moon. Un¬ 
fortunately, our CJS and EJS simulations did not record 
the collisional velocities between embryos, so the im¬ 
pact velocity distribution in these simulation sets is not 
known. However, we do have the collisional velocity data 
from the ANN simulations, which had 35 impacts be¬ 
tween Earth and Theia analogs. The cumulative distri¬ 


bution of impact velocities between Theia analogs and 
Ea rth analogs in ANN simulations is shown in Figure 
[TSl In this plot we see that the median impact velocity 
is just a few percent greater than the mutual escape ve¬ 
locity of the proto-Earth and Theia analogs with m > 0.1 
Mq. Moreover, the largest impact velocity seen in our 
simulations is 126% of the mutual escape veloci ty. For 
compari s on, m ost of the successful collisions in 


Cuk & 


Stewart (2012) require an impact velocity at least 150% 


ot the mutual escape speed. Thus, our ANN simulations 
suggest that such high velocity impacts are rare for Theia 
a nalogs with m > 0.1 M ^ 


Cuk & Stewart 


( 2012 [) also explore impact scenarios in- 


volving 'I'heia masses of 0.05 M 0 and 0.025 M 0 and find 
that such collisions can successfully produce an Earth- 
heavy composition for the Moon if the proto-Earth was 
spinning very rapidly before impact. Because of this find¬ 
ing, we also look at our collision statistics for last major 
mergers on Earth analogs that involve impacting bodies 
with masses below 0.1 M 0 . These are also shown in Fig¬ 
ure We see that smaller impactors do collide with 
the Earth at higher velocities, but the large majority are 
still below 150% of escape velocity, the minimum veloc- 


Cuk & Stewart 


(2012). Of the Theia 


ity preferred in 
analogs with 0.025 M 0 < m < O.iJ 5 M 0 , 17% strike the 
proto-Earth with a velocity that is more than 150% of 
the mutual escape velocity. For Theia analogs with 0.05 
M 0 < m < O.IM 0 , the percentage of high-velocity colli¬ 
sions falls to 7.9%. Because the CJS and EJS simulations 
involve collisions of bodies spanning a larger range of 
semimajor axes than the ANN set, we may expect larger 
fractions of energetic collisions in these two simulation 
sets. Unfortunately, t hese fractions are not k nown. How¬ 
ever, in order for the Cuk & Stewart (2012) mechanism 


to succeed, the Earth’s spin rate also has to be exception¬ 
ally high (a period of ~2.3 hours). Such a high spin rate 
is likely achieved with a very large merger event before 
the Moo n-fo rming impact, and we have already shown 
in Figure [TT] that massive impactors are rare. Thus, pro¬ 
duction of our moon with a high velocity impact may 
also be a relatively low probability event. 
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Fig. 19. — The distribution of | A^^O| predicted for Venus analogs 
if we assume that the initial distribution of varies linearly 

with heliocentric distance. The predicted values are shown for the 
CJS, EJS, and ANN simulations with the solid, dashed, and dotted 
CDFs, respectively. The shaded region marks 0.29%o< lA^^O] < 
0.32%o, the approximate observed value for Mars. 

3.5. Comparison to Venus Analogs 

With four different families of initial distribu¬ 

tions set by the accretion histories of Earth and Mars 
analogs, we can now predict the values that these 

would yield for Venus by examining the accretion his¬ 
tories of Venus analogs in our simulations. Because we 
need Earth and Mars analogs to define our initial A^^O 
distributions, we only use simulations that form analogs 
for Venus, Earth, and Mars. (This set of 70 simulations 
is slightly different from the last section since we drop 
the requirement for Theia analogs and replace it with 
Venus analogs.) As in the previous section we examine 
the outcomes of four different initial A^^O distributions 
for our embryos and planetesimals: a linear distribution, 
a step-function distribution, a random distribution, and 
two differe nt v alues for planetesimals and embryos. 

In Figure]^ we show the distributions of | A^^O| found 
for our Venus analogs if we assume the initial A^’^O 
among embryos and planetesimals varies linearly with 
their heliocentric distance. As before, we force the lin¬ 
ear distributions to yield |A^^O| ~ 0.32%o for the Mars 
analog in each simulation. When this is done, we find 
that the Venus analogs are unlikely to have similar oxy¬ 
gen isotope compositions to the Earth. In fact for each 
of our three simulation sets (CJS, EJS, and ANN), the 
median value of |A^^O| is only slightly more Earth-like 
than Mars. This is indicative of the unique feeding zones 
that Earth and Venus analogs have in a given system. 
Because Earth and Venus analogs have higher minimum 
masses than Mars analogs, it was suspected that the 
greater number of accretion events may cause the two 
planets to closely converge to the same average A^^O 
value. However, this is not the case. 

We now repeat this exercise, but instead the initial 
A^^O distribution is modeled with a step function with a 
given contrast and heliocentric position to yield 0.29%o < 
|A^^O| <0.35%o for the Mars analog of each system. For 
each chosen step function contrast, we calculate the me¬ 
dian value of |A^^O| for our Venus analogs in each en¬ 
semble (CJS, EJS, and ANN) if at least half of the sim¬ 
ulations yielded a Martian |A^^O| between 0.29%o and 
0.35%o. These median values are displayed in Figure [20^3. 




Fig. 20.— A: The fraction of Mars analogs that are able to 
satisfy |A^^O| = 0.32 di 0.03%o (compared to the system’s Earth 
analog) as a function of the chosen jump in |A^^O| in the initial 
|A^^O| step-function distribution (for all simulations containing 
Venus, Earth, and Mars analogs). CJS, EJS, and ANN simulation 
results are shown with black circles, blue diamonds, and red stars, 
respectively. B: The median |A^^O| value predicted for Venus 
analogs as a function of the chosen jump in |A^^O| imposed for 
the initial |A^^O| step-function distribution. CJS, EJS, and ANN 
simulation results are shown with black circles, blue diamonds, and 
red stars, respectively. 


As can be seen, the expected |A^^O| for Venus analogs is 
always at least O.I%o different from Earth (which by deh- 
nition is 0) and are often equal or greater than than Mar¬ 
tian |A^’^0| to which the initial distribution is tuned. It 
is also clear that the ANN simulations consistently yield 
the most Earth-like oxygen compositions for Venus. This 
is presumably because the ANN simulations more often 
yield Venus and Earth analogs that are much more mas¬ 
sive than the Martian analogs (Hansen 2009), and they 
more uniformly sample the narrow annulus of embryos 
and planetesimals compared to CJS and EJS simulations. 

Next we test the assumption that the initial A^^O of 
each embryo and planetesimal is uncorrelated with he¬ 
liocentric distance, and we draw them from a Gaussian 
distribution designed to typically yield a Martian |A^^O| 
near 0.32%o. This is repeated 10,000 times on each sim¬ 
ulation to build up robust statistics. Again, Figure [2lp 
shows that the |A^’^0| values for Venus analogs areWs- 
tinctly un-Earthlike. In fact, they once again often take 
values similar to Mars. 

Finally, we again explore the possibility that planetesi¬ 
mals are born with one fixed A^^O and embryos are born 
with another. As in section |3.3.4[ we set the two differ¬ 
ent values to yield A^^O = 0.32 %o for the IVIars analog 
in each simulation (relative to the Earth analog). Using 
this initial A^^O distribution we predict the A^^O val¬ 
ues for the Venus analogs of our simulations. These are 
shown in Figurej^ We see in this figure that this analy¬ 
sis predicts that Venus should be even more isotopically 
different from the Earth than Mars. The median |A^^O| 
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Fig. 21. — A: The distribution of |A^^0| for Mars analogs if we 
assume that the initial values of planetesimals are randomly 

drawn from a Gaussian distribution (for simulations that contain 
Earth, Mars, and Venus analogs). CJS, EJS and ANN simula¬ 
tions are shown with solid, dashed, and dotted CDFs, respectively. 
The shaded region marks 0.29%o< |A^’^0| < 0.32%o. Initial A^^O 
values of embryos and planetesimals were redrawn 10000 times to 
build our statistics. B: The distribution of | A^^O| values for Venus 
analogs for the same simulations as in panel A. The shaded region 
again marks 0.29%o< |A^^O| < 0.32%o. 



Fig. 22. — The distribution of | A^^O| predicted for Venus analogs 
if we assume two different values of A^^O for planetesimals and em¬ 
bryos. The predicted values are shown for the CJS, EJS, and ANN 
simulations with the solid, dashed, and dotted CDFs, respectively. 
The shaded region marks 0.29%o< |A^^O| < 0.32%o, the approxi¬ 
mate observed value for Mars. 

values for Venus are 0.73%c and 0.86%c for the CJS and 
EJS simulations, respectively. The reason for this is that 
the planetesimal mass fraction is almost always highest 
for the Earth, second highest for Mars, and third highest 
for Venus analogs in our simulations. Thus, Venus is al¬ 
most always predicted to be more dissimilar from Earth 
than Mars. 

Although Venus is nearly the same mass as Earth, we 
find that it is likely to have a significantly different 
value than the Earth. Thus, we find that Earth-size (and 


Venus-size) planets do not usually gravitate toward the 
same disk-averaged A^^O. This finding is unchanged for 
each of our four different initial A^^O distributions that 
we impose, as well as the three different sets of simulation 
initial conditions that we use. While this is by no means 
an exhaustive search of parameter space, it suggests that 
Earth and Venus should be expected to have distinct 
A^^O values. 

4. CONCLUSIONS 

We have undertaken a statistical study of the giant im¬ 
pact phase of terrestrial planet formation in order to un¬ 
derstand the relations between a planet’s final mass and 
orbit and its accretion feeding zone within the initial disk 
of embryos and planetesimals. We find that the width of 
a planet’s feeding zone does not correlate with its mass 
or semi-major axis, but is proportional to the width of 
the initial protoplanetary disk, with massive planets ex¬ 
hibiting remarkably little variance. This means that one 
should expect terrestrial planets of all sizes to sample the 
same parts of the disk, provided they form at the same 
distance from their star. In other words, super-Earths 
forming in the habitable zone of their star should have 
water mass fractions drawn from the same distribution 
as Earth-size planets. Moreover, although feeding zones 
are broad and often overlap in multi-planet systems, the 
constant feeding zone size dictates that planets forming 
farther from their star are more likely to accrete water- 
rich planetesimals from beyond the snow line. 

Eocusing on planets analogous to our own, we find that 
Venus analogs are typically built from material orbiting 
close to the Sun, and Mars analogs are constructed from 
more distant material, while Earth is a compromise be¬ 
tween these two accretion histories. However, these are 
only general trends. The wings of the feeding zones of 
each of these planets can span 1-3 AU, depending on 
the initial planetesimal disk and giant planet orbits we 
assume. In addition, there is a great deal of stochastic- 
ity in the planets’ feeding zones, which has a large effect 
on each planet’s inventory of distant, presumably water- 
rich, planetesimals. The differences in the planets’ outer 
tails of their feeding zones is largely ruled by random¬ 
ness, and it is not unusual for Venus to have a higher 
water mass fraction than Mars or Earth. 

Our simulations also allow us to study the compo¬ 
sition of Theia analogs relative to Venus, Earth, and 
Mars analogs, and this is a major goal of our work. 
We find that the feeding zones of Theia analogs vary 
wildly compared to the planets, and it cannot be assumed 
that Theia had a similar accretion history to Earth or 
any of the other planets. Consequently, it is unlikely 
that Theia’s oxygen isotope composition is similar to the 
Earth’s. Although we do not know the distribution of 
A^^O among our solar system’s initial planetesimals and 
embryos, we impose four different assumed A^^O dis¬ 
tributions designed to reproduce the observed A^’^0 of 
Mars. Regardless of which family of A^^O distributions 
we choose, we find that the probability of a Theia with 
|A^'^0| < 0.016%c is ~5% or less. Thus, a Theia with 
an oxygen isotope composition distinctly different from 
Earth’s is the expected outcome of terrestrial planet for¬ 
mation for any of the A^^O distributions we employ. A 
consequence of this is that our simulations typically pre¬ 
dict that the Moon cannot contain more than ~5-I5% 
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of Theia material, depending on our simulation set and 
our assumed initial distribution of A^^O. There may be 
more complex distributions that yield higher al¬ 

lowable Theia mass fractions for the Moon, but this is 
not the case for the relatively simple distributions that 
we employ. 

We focus on Theia’s oxygen composition since it has 
been the subject of so much past work, but the Moon 
is also known to be isotopically similar to t he Earth 
for tungsten, silicon, chromium, and titanium (j^ uboM 


et al. 2007 Armytage et al. 2012 Zhang et ar' 2012 ). 


Tresumably, the isotopes ot these other elements all had 
their own unique distributions in the protoplanetary disk 
that were different from oxygen. While we have shown 
that Earth and Theia can have feeding zones that by 
chance yield similar oxygen isotope compositions, the 
probability that these feeding zones would also simul¬ 
taneously yield similar isotope compositions for other el¬ 
ements is presumably much smaller. 

In addition, we examine the mass ratios between Theia 
analogs and the proto-Earth analogs they strike in our 
simulations. If Thei a ’s mas s is at least ~ 66 % that of the 
proto-Earth Canup (20121 demonstrates that the colli¬ 
sion of these two bodies can yield a moon with an Earth¬ 
like isotope signature even if the two bodies had dis¬ 
tinctly different initial compositions. However, we find 
that such collisions are extremely rare or absent in the 
final phase of terrestrial planet formation. With over 
100 collisions between Theia analogs and proto-Earth 
analogs, the largest mass ratio we observe is 0.48, and 
the large majority of collisions have mass ratios below 
0.4. 

Another scenario that can produce a terrestrial compo¬ 
sition for the Moon is a higher velocity impact between 


the proto-E arth and Theia (Cuk fc Stewart|2012 


Reufer 


et al. 2012 1 . In one set of our simulations we examine 
the distributions of impact velocities between Earth and 
Theia analogs, and we find that for Theia analogs with 
TO > 0.1 M 0 the typical impact velocity is just a few per¬ 
cent greater than the mutual escape velocity. Meanwhile, 
most high velocity collisions that successfully reproduce 
the moon require an impact speed that is at least ^50% 
greater than the escape speed. If we consider Theia ana¬ 
log masses down to 0.025 M 0 there are a significant num¬ 


ber of collisions above 150% of the escape speed, but this 
is still not the most common outcome. Furthermore, the 
proto-Earth must presumably be spun up by a massive 
collision before the moon-forming impact, and we find 
these massive collisions to be rare. This suggests that a 
high velocity moon-forming collision may also not be a 
high probability event. 

Our work demonstrates that although the giant impact 
hypothesis is the most favored scenario for the origin 
of the moon, the isotopic similarity between the Earth 
and the Moon conflicts with predictions of terrestrial 
planet formation simulations using the simple A^^O dis¬ 
tributions employed here. One potential solution to this 
is that the outer Earth and moon-forming disk under¬ 
went substantial mixing immediat ely after Thei a’s im¬ 
pact (i Pahlevan fc Stevenson||2007|) . However, as |Canup| 
(20131 points out, m this scenario one would expect 
larger isotopic differences between the Earth and the 
Moon to exist for more refractory elements. Such a trend 
has not been shown at present. As a result, all explana¬ 
tions for the Moon and Earth’s isotopic similarities seem 
to rely on low probability events, and perhaps our Moon 
is a statistical outlier in the context of giant impacts. 

Finally, using the oxygen isotopic distributions we have 
generated to study Theia analogs, we also make predic¬ 
tions for the composition of Venus analogs in our simu¬ 
lations. At present there are no known samples of Venus 
material. However, the oxygen isotope distributions we 
invoke to explain the isotopic signature of Mars predict 
that Venus will also have an isotope composition dis¬ 
tinct from Earth, regardless of which initial isotope dis¬ 
tribution we choose. The |A^^O| values we predict for 
Venus are comparable to or greater than those observed 
for Mars, depending on our assumed initial A^^O distri¬ 
bution. 


5. ACKNOWLEDGEMENTS 

This work was partially funded by NSF Grant AST- 
1109776. We are very grateful to John Chambers, Sean 
Raymond, and Erik Petersen for helpful discussions while 
preparing our work for publication. In addition, com¬ 
ments and suggestions from the two anonymous referees 
greatly improved this paper. 


REFERENCES 


Abbot, D. S., Cowan, N. B., & Ciesla, F. J. 2012, ApJ, 756, 178 
Armytage, R. M. G., Georg, R. B., Williams, H. M., & Halliday, 
A. N. 2012, GeCoA, 77, 504 

Burbine, T. H., & O’Brien, K. M. 2004, Meteoritics and 
Planetary Science, 39, 667 

Cameron, A. G. W., & Ward, W. R. 1976, in Lunar and 

Planetary Inst. Technical Report, Vol. 7, Lunar and Planetary 
Science Conference, 120 
Canup, R. 2013, Nature, 504, 27 
Canup, R. M. 2004, Icarus, 168, 433 
—. 2012, Science, 338, 1052 

Canup, R. M., & Asphaug, E. 2001, Nature, 412, 708 
Chambers, J. E. 1999, MNRAS, 304, 793 
—. 2001, Icarus, 152, 205 

Clayton, R. N., & Mayeda, T. K. 1996, GeCoA, 60, 1999 
Cowan, N. B., & Abbot, D. S. 2014, ApJ, 781, 27 
Cuk, M., h Stewart, S. T. 2012, Science, 338, 1047 
Dauphas, N., Burkhardt, C., Warren, P., & Teng, F.-Z. 2014, 
Royal Society of London Philosophical Transactions Series A, 
372, 244 


Fischer, R. A., Ciesla, F. J. 2014, Earth and Planetary Science 
Letters, 392, 28 

Franchi, L A., Wright, L P., Sexton, A. S., k, Pillinger, C. T. 

1999, Meteoritics and Planetary Science, 34, 657 
Hansen, B. M. S. 2009, ApJ, 703, 1131 
Hartmann, W. K., k Davis, D. R. 1975, Icarus, 24, 504 
Hayashi, C. 1981, Progress of Theoretical Physics Supplement, 

70, 35 

Herwartz, D., Pack, A., Friedrichs, B., k Bischoff, A. 2014, 
Science, 344, 1146 

Izidoro, A., de Souza Torres, K., Winter, O. C., k Haghighipour, 
N. 2013, ApJ, 767, 54 

Izidoro, A., Haghighipour, N., Winter, O. C., k Tsuchida, M. 
2014, ApJ, 782, 31 

Jacobson, S. A., k Morbidelli, A. 2014, Royal Society of London 
Philosophical Transactions Series A, 372, 174 
Jacobson, S. A., Morbidelli, A., Raymond, S. N., O’Brien, D. P., 
Walsh, K. J., k Ruble, D. C. 2014, Nature, 508, 84 
Javoy, M., et al. 2010, Earth and Planetary Science Letters, 293, 
259 





































16 


Kasting, J. F., Whitmire, D. P., & Reynolds, R. T. 1993, Icarus, 
101 , 108 

Morishima, R., Schmidt, M. W., Stadel, J., & Moore, B. 2008, 
ApJ, 685, 1247 

O’Brien, D. P., Morbidelli, A., & Levison, H. F. 2006, Icarus, 184, 
39 

Ozima, M., Podosek, F. A., Higuchi, T., Yin, Q.-Z., & Yamada, 
A. 2007, Icarus, 186, 562 

Pahlevan, K., & Stevenson, D. J. 2007, Earth and Planetary 
Science Letters, 262, 438 

Quarles, B., & Lissauer, J. J. 2014, ArXiv e-prints 
Quintana, E. V., & Lissauer, J. J. 2014, ApJ, 786, 33 
Raymond, S. N., O’Brien, D. R, Morbidelli, A., & Kaib, N. A. 
2009, Icarus, 203, 644 

Raymond, S. N., Quinn, T., & Lunine, J. 1. 2004, Icarus, 168, 1 
—. 2005, ApJ, 632, 670 
—. 2006, Icarus, 183, 265 


Reufer, A., Meier, M. M. M., Benz, W., & Wieler, R. 2012, 
Icarus, 221, 296 

Touboul, M., Kleine, T., Bourdon, B., Palme, H., & Wieler, R. 
2007, Nature, 450, 1206 

Walsh, K. J., Morbidelli, A., Raymond, S. N., O’Brien, D. P., & 
Mandell, A. M. 2011, Nature, 475, 206 
Wiechert, U., Halliday, A. N., Lee, D.-C., Snyder, G. A., Taylor, 
L. A., & Rumble, D. 2001, Science, 294, 345 
Zhang, J., Dauphas, N., Davis, A. M., Leya, L, & Fedkin, A. 
2012, Nature Geoscience, 5, 251 

Ziegler, K., Ghambers, J. E., & Young, E. D. 2006, in Lunar and 
Planetary Inst. Technical Report, Vol. 37, 37th Annual Lunar 
and Planetary Science Conference, ed. S. Mackwell & 

E. Stansbery, 1894 



